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ABSTRACT 

We investigate generation processes of magnetic fields around first stars. Since 
the first stars are expected to form anisotropic ionization fronts in the surround- 
ing clumpy media, magnetic fields are generated by effects of radiation force as 
well as the Biermann battery effect. We have calculated the amplitude of mag- 
netic field generated by the effects of radiation force around the first stars in 
the preceding paper, in which the Biermann battery effects are not taken into 
account. In this paper, we calculate the generation of magnetic fields by the 
Biermann battery effect as well as the effects of radiation force, utilizing the 
radiation hydrodynamics simulations. As a result, we find that the generated 
magnetic field strengths are ~ 10 _19 G — 10~ 17 G at ~ lOOpc-lkpc scale mainly by 
the Biermann battery, which is an order of magnitude larger than the results of 
our previous study. We also find that this result is insensitive to various physical 
parameters including the mass of the source star, distance between the source 
and the dense clump, unless we take unlikely values of these parameters. 

Subject headings: early universe — HE regions — radiative transfer — magnetic 
fields 



1. Introduction 

According to the theoreti cal studies in the last decade, first stars are expected to be 



very massive (> lOOM^) (e . g.. iBromm. Coppi fc Larsonll2002l ; iNakamura fc Umemural 12001 



Abel. Bryan. Norman 2002 ; Yoshida 2006h . Recent studies which properly address the ac- 



cretion phase of first star formation also revealed that the primary star formed in the center 
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of the mini-halo is not very massive but still massive > 10 M^, although signi fic ant fraction 



of fir st stars are less massive (< lM0) (jStacy et al.ll2010l ; IClark et al.ll2011al Jbl; iGreif et al. 



201ll ). In any case, star formation episodes in the very early universe are different from that 
in local galaxies. One of the reasons of this difference is that the primordial gas clouds that 
host the first stars lack heavy elements, though they are most efficient coolants in interstellar 
clouds at T < 1000K. Because of the lack these coolants, the temperature of the primordial 
gas cloud is kept around ~1000K during its collapse for % < 10 12 cm~ 3 , which is much 
hotter than the local interstellar molecular gas clouds. Consequently, the Jeans mass of the 
collapsing primordial gas is much lar ger than that of the interstellar gas, which leads to the 
formation of very massive starsfe.g.. IOmukaill2000l ). 



Another important difference between the star formation sites in the early universe and 
local molecular clouds is the strengths of magnetic fields. Typical field strength B ~ 10/iG in 
the local molecular gas results in the formation of jets from protostars and regulate the gravi- 
tational collapse of cloud cores. The effects of magnetic field on the star formation in the early 
universe have been studied from theoretical aspects. First of all, the coupling of the magnetic 
field with t he primordial gas was studied by a detailed chemical reaction mode KjMaki fc Susa 
2004 . 120071 ) . They found that magnetic field is basically fr ozen-in the primordial gas dur ing 
its collapse, differently from the local interstellar gas(e.g., Nakano fc Umebayashi 1986al lbh . 
Under the assumption of the flux freezing condition, the dynamical importance of the mag- 
netic field is also investigated by several authors. In case the magnetic field is stronger than 
~ 10~ 9 G at tt-h = 10 3 cm -3 , field strength is amplified to ~ 10 3 G at wh = 10 21 cm -3 which is 
eno ugh to launch the bipolar o ut flows, and to suppress th e fragmentation of the accretion 
disk OMachida et al.l 120061 . 120081 ) . iTan fc Blackmanl ( 120041 ) estimated the condition for the 
magnetorotational instability (MRI) to be activated in the accretion disk around the proto- 
star, by comparing the Ohmic dissipation time scale with the growth time scale of MRI. They 
found that the condition is B > 10~ 2 G at hh = 10 15 cm~ 3 which corresponds to B > 10~ 10 G 
at wh = 10 3 cm -3 . We also remark that the turbulent motion powered by the accreting gas 
can amplify the initial field strength faster than the simple flux freezing, although the effects 
are still under debate. Thus, the magnetic field could be of importance at the final phase of 
the star formation process in primordial gas clouds, if B > 10~ 10 — 10 _9 G at nn = 10 3 cm -3 , 
i.e. at the initial phase of the collapse of primordial gas in the mini-halos. 

In addition, recent theoretical studies suggested that the heating by the ambipolar dif- 
fusion process in star-forming gas clouds could change the thermal evolution o f the prestellar 
core in case the fie ld stre ngth is as strong as 10 _10 G at IGM comoving densities (ISchleicher et al 



2009 



Sethi et al. 



2008). This process also might leads to the formation of massive black 



holes ( Sethi et al. 201ol ) sin ce such heating can shu t down the H2 cooling and open the path 



of the atomic coolingfe.g.. lOmukai & Yoshiil 120031 ) . 
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In any case, it is important to determine the magnetic field strengths in star forming 
gas clouds in the early universe, in order to quantify the effects of magnetic fields on the 
primordial star formation. In spite of such potential importance, initial seed magnetic field 
strengths are still unknown observationally. Only the observations on the distortion of the 
cosm ic microwave background spectrum(e.g., iBarrow et al.lll997l ; ISeshadri fc Subramanian 



20091 ). and the measu rements of the Faraday rotati on in the polarized radio emission from 



distant quasars (e.g., iBlasi et al.lll999l ; IValled 12004 ) imply rather mild upper limits on the 
field strengths in IGM, at the level of B ~ 10~ 9 G. On the other hand, various theoretical 
studies predicted that it is as small as < 10~ 18 G at IGM densi ties. For instanc e, there are 
models generating th e magnetic field by t he Biermann b attery (jBiermannlll950l ) during the 
structure f ormat ion (IKulsrud et al.lll997l ; Ku et al.l 120081 ). Recent numerical simulation by 



Sur et al.l (120101 ) suggests that strong magnetic field emerges during collapse of turbulent 



prestellar cores of primordial gas due to the Biermann battery and the turbulent dynamo 
action. T here is also a number of models that the fields are generated just after the big 
bang(e.g., Ilurnerfe Widrowlll988l ; llchiki et al.l 120061 ) . 

It is also s ugges ted that reionization of the universe inevitably generates magnetic fields. 



Gnedin et al.l ( 120001 ) have shown that the considerable Biermann battery term arises at the 
ionization fronts in their cosmological simulations. They predict ~ 10 _18 G a t Sp/p ~ 10 3 . 



It als o is suggested that in the neigh borhood of luminous sources like QSOs fLanger et al. 



20031 ) or first stars ( lAndo et al.l ( 120101 ); here after ADS10) magnetic field could be generated 



through the momentum transfer process from ionizing photons to electrons, since V x E ^ 
is satisfied at the borders between the shadowed regions and ionized regions. ADS10 
predicted B ~ 10~ 19 G at IGM densities at z = 20, however, they did not take into account 
the Biermann battery effect, since they assume that the gas is isothermal and static. 

In this paper, we extend our previous study (ADS10) to investigate the generation pro- 
cess of magnetic fields due to the ionizing radiation from first stars with more precision. We 
take into consideration not only the effects of radiation force but also the Biermann battery 
mechanism, utilizing the two-dimensional radiation hydrodynamics simulations. Then we 
discuss whether the magnetic field strength obtained in our study could be important for 
subsequent star formation process. In section 2, we describe the basic equations and the 
setup of our model. We show the results of our calculations in section 3. Sections 4 and 5 
are devoted to the discussions and summary. 
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2. Basic equations &c Model 

2.1. Equation of magnetic field generation 

According to ADS 10, the equation of magnetic field generation is given as 
OB c c 

— = V x (v x B) -Vn c x V Pc - -V x / rad (1) 

at erig e 

where .B, v and e denote the magnetic flux density, fluid velocity, and the elementary charge, 
respectively. p e and n e represent the pressure and the number density of electrons. / ra d is 
the radiation force acting on an electron. The first term on the right-hand side is the 
advec tion term of mag netic flux, whereas the second term describes the Biermann battery 



term (IBiermannl 119501 ) . which was not included in our previous study (ADS 10). The third 
term is the radiation term which represents the momentum transfer from photons to gas 
particles. Remark that the dissipation term due to the resistivity of the gas is omitted, since 
it is negligible in comparison with the other terms (ADS10). 



2.2. Radiation force and photoionization rate 

The radiation force on an electron, / ra d involves two processes. The first one is the 
contribution by the Thomson scattering, / ra d,T- /rad,T is given by 

»1>L 



/rad,T = — / F 0u dis-\ / F 0u exp[-T UL a(is)]du. 

c JO c Ju L 



(2) 



where o"t denotes the cross section of the Thomson scattering, F 0u is the unabsorbed energy 
flux density, ul is the Lyman-limit frequency, r Vt denotes the optical depth at the Lyman 
limit regarding the photoionization, and a (is) is the frequency dependence of photoionization 
cross section which is normalized at the Lyman limit frequency. 

Another source of the radiation force is the momentum transfer from photons to electrons 
through the photoionization process. The force expressed as / ra d,i, is 

/rad,l = r — / 0v L a 0) F 0u exp [-r UL a (is)} dis, (3) 
2 cn c J UL ' - 

where a VL is the photoionization cross section at the Lyman limit and nni represents the 
number density of neutral hydrogen atoms. 
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In order to assess the electron number density, we also solve the following photoionization 
rate equation for electrons: 

Qy e r 

-t^t + {v ■ V)y c = a B y e y p n u + k coll y e y m n R , (4) 

at n H 

where «b and fc co u denote the case B recombination rate and collisional ionization rate per 
unit volume, respectively. y e , y p and ym is the number fraction of electrons, protons and 
neutral hydrogen atoms, respectively, nn is the number density of hydrogen nuclei and 
v denotes the velocity of the fluid. The photoionization rate per unit volume, T, is also 
obtained by the formal solution of the radiation transfer equation: 

a UL -j^exp \-T VL a (y)] du. (5) 



2.3. Hydrodynamics with heating/cooling 

We solve the ordinary set of hydrodynamics equations: 

: -pV-V, (6) 



and the equation of state 



Dp 
Dt 

= _Ivp-V0 DM , (7) 
Dt p 

P^ t = G-L-pV-v, (8) 



V = ( T - l)pe, (9) 



where p, p, and e are the density, pressure, and the specific energy of the fluid, respectively. 
7 denotes the specific heat ratio. G and L are the radiative heating rate and cooling rate 
per unit volume, respectively. 0dm is the gravitational potential of the dark matter halo(see 
12. 4ft . The feedback from magnetic fields to the fluid is neglected in equation^, since we 
consider the generation of very weak magnetic fields. We also omit the self gravitational 
force of gas, which is unimportant as long as we consider the gas with nn < 10 3 cm~ 3 . 

In order to perform hydrodyna mics simulations, we use the Cubic Interpolated Profile 



(CIP) method OYabe fc Aokilll99ll ). The CIP scheme basically tries to solve not only the 
advection of physical quantities but also the derivatives of the quantities. Using this scheme, 
we can capture spatially sharp profiles of fluids, that are always expected in the problems 
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including the propagation of ionization fronts. In this paper, the CIP scheme is applied to 
the advection terms of equations (TjJ and (jSJ)-©- 

Using the formal solution of radiation transfer equations, the radiative heating rate G 
is assessed as 



G = n m 



OO 771 



hv 



h{y — ul) o VL a iy) exp [—r UL a (y)] dv. 



(10) 



The cooling rate, L, is given as 

L = L co nn e n m + L vec n e n p + L exc n e n m + L s n e n p , 



where L co \\, L rec , L exc and Lg are the cooling coefficients regarding the collisional-ionization, 
the recombination, the collisional e xcitation and the free- f ree em ission. These cooling rates 
are taken from the compilations in iFukugita fc Kawasaki! (11994J ). 



2.4. Setup 



We consider a mini-halo in the IGM at redshift z ~ 20 exposed to an intense radiation 
flux from a nearby first star. We assume that the dark matter density profile of the mini-halo 
is described by the NFW profile (jNavarro et al.lll997l ) 



Pdm \r) 



(r/r s ) (1 + r/r s 



\2 " 



121 



where r is the radial destance from the center of halo. p s and r s are a characteristic 
density and radiu s, which are determined by a halo mass Mh a i collapsing at redshift z 
(IPrada et al.ll201ll ). We assume the gas density profile is a core-halo structure with p oc r~ 2 
envelope. The core radius r c is determined by the core density and the total gas mass 
(M gas = (Qb/^m) Mhaio)- We study various cases of core densities (n cm -3 ), the halo mass 
M ha io and distances between the source star and halo center (Dkpc). We assume stationarity 
of the mini-halo with respect to the source star. This assumption is based upon the fact 
that the the change of distance between the star and the halo due to the relative motion is 
smaller than D, if we consider the cosmic expansion. We remark that in case the neighboring 
overdense region and the source star is contained in a same halo of > 10 7 M Q , change of the 
distance due to the velocity dispertion of the halo could have significant effect especially for 
low mass source stars. 
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Initially, the ambient gas is assumed to be neutral when the source star is turned on. The 
initial number density of the ambient gas in the intergalactic space is uiqm = 10~ 2 cm~ 3 . We 
assume that the initial temperature is 500K. We perform two-dimensional simulations in the 
cylindrical coordinates assuming axial symmetry. As shown in Figure dj the computational 
domain is 100 pc x 200 pc in R — z plane. We consider source stars of various masses, 
M* = 500,300,120,60, 25 and 9M . The luminosit i es, ef fective temperatures and the ages 
of these stars are taken from the table of ISchaererl (120021 ). The incident radiation from the 
source star is assumed to be perpendicular to the left edge of the computational domain. 

We employ four models (A-D) of different Mh a io,^o an d D listed in Table [fl these are 
plausible values of the minihalos in standard ACDM cosmology (see section [372]) . The number 
of grids we use in these simulations is basically 250 x 500 which is confirmed to be enough 
to obtain physical results by the convergence study (see section | 



3. Results 

3.1. Typical results 

First, we consider a first star of 500 M & (t age = 2 x 10 6 yr). Figure [2] shows the results 
for the model A. Two columns correspond to the snapshots at O.lMyr and 2Myr. The top 
row shows the color contour of the mass density of the gas. The gas density distribution is 
isotropic at O.lMyr, while it is highly disturbed by the radiation flux from the left at 2Myr. 
We also can find a shock front at 2Myr inside the core of the overdense region. 

The second and third rows show the maps of the temperature and the electron density. 
Clearly, an ionization front is generated by the flux from the left, and it propagates into 
the core of the cloud. It is also worth noting that the pattern of temperature and electron 
density distribution are similar, but not identical with each other. Slight difference between 
the two contour maps directly leads to the nontrivial Biermann battery term. 

The bottom row shows the distribution of the generated magnetic field strength. As 
expected, we obtain strong magnetic fields in the neighborhood of the ionization front. In 
addition, strong magnetic fields is produced at the center of core where the density is the 
highest. The peak field strength is as strong as 6 x 10~ 18 G in this case. 

In Figure [3j we show the time evolution of the peak magnetic field strength in the 
simulated region (red curve). We also plot the peak magnetic field generated by the Biermann 
battery term (blue) as well as the one by the radiation processes (green). The field strength 
grows almost linearly in this model, and the final strength is as large as ~ 6 x 10~ 18 G. 
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We also find that the radiation process is less important than the Biermann battery effect. 
However, it is still noteworthy that the difference between the two contributions is only a 
factor of ~ 10, although the nature of these processes are very different from each other. 

3.2. Dependence on n , D, 

We also investigate the other models of n and D listed in Table [TJ The snapshots at 
2Myr of these models (B,C and D from top to bottom) are shown in Figure HI The left 
column shows the magnetic field strength, whereas the right column illustrates the number 
density of electrons. In the model B, the generated magnetic field is of the order of ~ 10~ 19 G 
(top left panel), which is smaller than that in the model A. Since the ionization front is not 
trapped by the dense core in the model B (top right), the source terms of equation flTJ 
become very small as soon as the ionization front passes through the core. As a result, the 
magnetic field does not have enough time to grow. In the models C and D, the core is 10 
times more distant from the source star than the models A and B. Therefore, the ionization 
front is trapped at lower density regions, and becomes less sharp than that in model A. 
Consequently, the generated magnetic field is smaller than that in the model A (middle and 
bottom row). 

We also plot the time evolution of the peak field strengths of the models A-D in Figure 
[5j The models A, C and D mostly increase monotonically, whereas the model B has a clear 
plateau/decline. Such different behavior also comes from the fact that the ionization front 
immediately passes through the core in the model B. In such case, the time for the magnetic 
field to grow is not enough as stated above. In addition, the fluids that host the generated 
magnetic fields expand due to the thermal pressure of photoionized gas. Such expansion 
results in the slight decline of the magnetic field strength. 

The dependence on the source stellar mass M* is also studied. We employ six models 
of M* = 500,300,120,60,25 and 9M , while the other parameters are same as the model A. 
Figure [6] shows the peak magnetic field strength as a function of M*. The peak field strength 
of each run is evaluated when it gets to the time for the death of the source star. The peak 
field strength basically increases as the stellar mass becomes more massive. However, the 
difference between the field strength of M* = 9M and 500M Q is a factor of 9, which is 
a small difference for such a mass difference. The reason of this behavior comes from the 
fact that the more massive the stars are, the more ionizing photons they emit, but also the 
shorter lifetime they have. These two competing effects cancel with each other. 

In any case, the generated magnetic field strengths stay around 10~ 18 — 10~ 17 G, if we 



-9 - 



consider the stellar mass of 9M & - 5OOM , which range is wide enough for first stars. 



The parameters n and D we employ in this paper are reasonable values. The baryonic 
density of the virialized halo at z = 20 is ~ 1cm -3 , and the sizes of the first halos are 
< lOOpc. In ad dition, the radii of c osmological HII regions in the numerical simulations are 
a few kpc (e.g.. lYoshida et al.ll2007l ). Thus, our models A-D are reasonable for the standard 
cosmological model. 



To summarize the results of this section (13. 2p . the magnetic field strength generated 
around a first star is 10~ 19 G < B < 10~ 17 G for M* = 5OOM , and a factor of a few smaller 
for less massive stars. 



3.3. Coherence length 

In addition, we also remark that the coherence length of the magnetic field. Due to the 
limited computational resource, we employ rather small box size (~ 200pc). The resultant 
coherence length of the magnetic field clearly exceeds the box size in most of the cases. In 
addition, if the ionization structure is more or less similar to that of our previous results in 
ADS10, the coherence length will be as large as ~ lkpc, which is the size of the shadow. 



3.4. Convergence check 



Since the equation of the magnetic field generation (JTJ includes spatial derivatives in its 
source terms, we have to pay attention to the effects of the cell size on our numerical results. 
We check the numerical convergence of the magnetic field strengths for the model A. We 
perform runs with N R x N z = 63 x 125, 125 x 250, 250 x 500 (canonical), 400 x 800, 500 x 1000, 
where Nr,Nz is the number of grids in i?-axis and z-axis, respectively. In Figure the 
magnetic field probability distribution functions (PDFs) of these runs are plotted. The axis of 
abscissas is the absolute value of the magnetic field strength, while the vertical axis represents 
the fraction of grid cells that fall in the range of [B, B + AB], where AB = 5 x 10 _19 G. The 
peak field strength decline as the number of grids increases. However, the peak field strength 
converges ~ 6 x 10~ 18 G and PDFs show similar distributions for Nr x Nz > 250 x 500. 
Therefore, the numerical results of the simulations converge very well at Nr x Nz = 250 x 500, 
which we use in other runs. 
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Discussions 



In this paper, we investigate the magnetic field generated by first stars. As a result, the 
maximal magnetic field strength is < 10~ 17 G, mainly generated by the Biermann battery 
mechanism. 

In fact, the order of magnitude of the magnetic field generated by the Biermann battery 
could be assessed as 

B ~^G£) (£) Sin ^ age 



,0 x lO^G ( J^-) ( S ^f) 
\2MjtJ V 0-1 J 



Ar TV T \ 

Ipc J \Wk) 



(13) 



where Ar denotes the typical length of n e and p e change significantly, 6 is the angle between 
Vn e and VT. The resultant value of above equation is roughly consistent with the results 
of our numerical simulations. 

We find the relative importance of the Biermann battery effect versus the radiative 
processes in this paper. However, this is only relevant for present setup, because the two 
effects depend differently on various parameters. In particular, the Biermann battery effect 
do not depend on the flux of the source star directly (see eq JT3|) . whereas the radiation force 
is proportional to the flux. This means if consider the magnetic field generation process in 
the very neighbor of the source objects, such as the accretion disks of the protostars/black 
holes, radiative processes could play central roles in the generation of magnetic field. We 
will study this issue in the near future. 



The field strength obtained in this paper is similar to the results of Ku et al.l ( 120081 ). 
in which they investigated magnetic fields in collapsing mini-halos/prestellar cores. If we 
assume that the magnetic field generated around first stars are brought into another prestellar 
core, and evolve during the collapse of the primordial gas in a similar way to Ku et al. 
( 120081 ). the magnetic field will be amplified up to ~ 10~ 13 G at 10 3 cm -3 . However, this 
magnetic field will not af fect subsequent star formation, since the magnetic field strength 
required for jet formation jMachida et al. 2006 ) and MRI activation( Tan fc Blackman 2004 ) 
is 10~ 10 — 10~ 9 G at 10 3 cm -3 , if we assume simple flux freezing condition. On the other 
hand, recent studies suggest t hat weak seed magnet i c field is ampli fied by the turbulenc e 
during the first star formation (ISchleicher et al.ll2010l ; ISur et al.ll2010l ). In ISur et al.l ( 120101 ) . 
weak seed magnetic fields are exponentially amplified by small-scale dynamo action if they 
employ sufficient numerical resolutions. In this case, the magnetic fields generated by the 
first stars in this paper might be amplified and affect subsequent star formation. However, 
these studies assume a priori given turbulence and initial magnetic field. To try to settle 
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this issue, we need cosmological MHD simulations with very high resolution. 

5. Summary 

In summary, we have investigated the magnetic field generation process by the radiative 
feedback of first stars, including the effects of radiation force and the Biermann battery. As 
a result, we found 10~ 19 G < B < 10~ 17 G on the boundary of the shadowed region, if we take 
reasonable parameters expected from standard theory of cosmological structure formation. 
The resultant field strength with a simple assumption of flux freezing suggests that the such 
magnetic field is unimportant for the star formation process. However, it could be important 
if the magnetic field is amplified by turbulent motions of the star forming gas could. 

We appreciate the anonymous referee for helpful comments. We also thank N.Tominaga 
and M.Ando for fruitful discussions. This work was supported by Ministry of Education, 
Science, Sports and Culture, Grant-in-Aid for Scientific Research (C), 22540295. 
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Model 


A 


B 


C 


D 


no [cm -3 ] 
L>[kpc] 
Magnetic field [G] 


5 x 10 6 
10 
0.2 

5.7 x 10~ 18 


1 x 10 5 
1 

0.2 

1.5 x HT 19 


5 x 10 6 
10 
2 

8.4 x 10~ 19 


1 x 10 5 
1 
2 

2.1 x 10~ 18 



Table 1: The peak magnetic field strengths of four models of M* = 500M ( 
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Fig. 1. — Schematic view of the computational domain. The shadowed region is formed 
behind the over dense region , while the other region is exposed to radiation field. 
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Fig. 2. — Two snapshots (left:0.1Myr, right:2Myr) for the model A are shown. Four rows 
correspond to the mass density of gas [gem -3 ] (top), the gas temperature [K], the number 
density of electron [cm -3 ], and the magnetic field strength [G], respectively. Black solid lines 
represent the position of the ionization fronts. 
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Fig. 3. — Time evolution of the peak magnetic field for the model A (red). Other two curves 
correspond to the field generated by Biermann battery (blue) and radiation pressure (green). 
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Fig. 4. — Snapshots at 2Myr for the models B,C and D are shown. Left column : magnetic 
filed strength, B[G] . Right column: electron density, n e [cm~ 3 ]. 
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Fig. 5. — Time evolution of the peak magnetic field for the models A-D. 
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Fig. 6. — The maximal magnetic field strength for D = 200pc, uq = 10cm -3 , with various 
M*. 
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Fig. 7. — Magnetic field probability distribution functions are plotted for various Nz- 



